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Abstract 

There has been much interest in possible violations of Lorentz invariance, partic- 
ularly motivated by quantum gravity theories. It has been suggested that a small 
amount of Lorentz invariance violation (LIV) could turn off photomeson interactions 
of ultrahigh energy cosmic rays (UHECRs) with photons of the cosmic background 
radiation and thereby eliminate the resulting sharp steepening in the spectrum of 
the highest energy CRs predicted by Greisen Zatsepin and Kuzmin (GZK). Recent 
measurements of the UHECR spectrum reported by the HiRes and Auger collabo- 
rations, however, indicate the presence of the GZK effect. We present the results of 
a detailed calculation of the modification of the UHECR spectrum caused by LIV 
using the formalism of Coleman and Glashow. We then compare these results with 
the experimental UHECR data from Auger and HiRes. Based on these data, we 
find a best fit amount of LIV of x 10 -23 , consistent with an upper limit of 

6 x 10 -23 . This possible amount of LIV can lead to a recovery of the cosmic ray 
spectrum at higher energies than presently observed. Such an LIV recovery effect 
can be tested observationally using future detectors. 
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1 Introduction 


Because of their extreme energy and isotropic distribution, it is believed that 
UHECRs are extragalactic in origin. After the discovery of the cosmic back- 
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ground radiation (CBR), Greisen [1] and Zatsepin and Kuzmin [2] pointed 
out that photomeson interactions should deplete the flux of cosmic rays with 
energies above ~ 50 EeV. One of us [3], using data on the energy depen- 
dence of the photomeson production cross section, then made a quantitative 
calculation of this “GZK effect” deriving the mean photomeson energy loss 
attenuation length for protons as a function of proton energy. These results 
indicated that the attenuation length of a proton with an energy greater than 
100 EeV is less than 100 Mpc, which is much less than the visible radius of 
the universe. Thus, what is sometimes referred to as the GZK “cutoff” is not 
a true cutoff, but a suppression of the ultrahigh energy cosmic ray flux aris- 
ing from a limitation of the proton propagation length through the cosmic 
background radiation owing to energy losses. 

From time to time there have been reports in the literature of the detection 
of giant air shower events from primaries with energies above the GZK sup- 
pression energy (trans-GZK events) ( e.g ., Refs. [4] - [6]). Such events have 
stimulated suggestions that a violation of Lorentz invariance or a modification 
of the Lorentz transformation relations at ultrahigh energies could result in 
a nullification of the GZK effect [7], [8]. Most significantly, the AGASA group 
reported 11 events above the GZK suppression energy [6], increasing the in- 
terest in the possibility of such new physics [9]. See Ref. [10] for a recent review 
of this topic. 

However, a reanalysis of the AGASA data (unpublished) has resulted in cut- 
ting their originally reported number of trans-GZK events by half. More im- 
portantly, the HiRes [11] and Auger groups [12], with larger exposures, have 
very recently claimed to have found a GZK suppression effect. Motivated by 
these new results, we have undertaken new detailed calculations of the effect 
of a very small amount of Lorentz invariance violation (LIV) on the spectrum 
of LIHECRs at Earth. We present our results here and compare them with the 
HiRes and Auger data separately. 


2 Violating Lorentz Invariance 


With the idea of spontaneous symmetry breaking in particle physics came 
the suggestion that Lorentz invariance (LI) might be weakly broken at high 
energies. Some of the more recent motivation for LIV has been in relating it to 
possible Planck scale phenomena that could lead to astrophysically observable 
consequences [13]. The idea that Planck scale physics may lead to a natural 
abrogation of the GZK effect has been of particular interest, since this would 
lead to a direct observational test. Significant fluxes of UHECRs at trans-GZK 
energies, could be the result of a very small amount of LIV [14] - [18]. Such 
a test would have important implications for some quantum gravity and large 
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extra dimension models, since those models may predict very small amount 
of LIV. 

Although no true quantum theory of gravity exists, it is natural to tie LIV to 
various quantum gravity models. A few examples of such work can be found 
in Refs. [15] - [20]. For more references, we refer to an excellent review by 
Mattingly [21]. A data table of constraints on LIV and CPT violation param- 
eters within the framework of the “Standard Model Extension” model [22] has 
recently been given by Kostelecky and Russell [23]. 

In this paper, we reinvestigate the observational implications of the possible 
effect of a very small amount of LIV, viz., that cosmic rays could indeed reach 
us after originating at distances greater than 100 Mpc without undergoing 
large energy losses from photomeson interactions. We considered this topic 
before in a more simplistic manner [24] when there was a clear discrepancy 
between the AGASA group data [6] and the earlier HiRes data. However, as 
discussed above, the observational situation has changed and now requires a 
more detailed approach. We therefore undertook a detailed calculation of the 
modification of the LIHECR spectrum caused by LIV using the formalism of 
Coleman and Glashow and the kinematical approach originally given by Alfaro 
and Palma [18] in the context of the Loop Quantum Gravity model [25], [26]. 
(See also Ref. [27].) Then, by comparing our results with the observational 
UHECR data we can place a quantitative limit on the amount of LIV. We also 
discuss how a small amount of LIV that is consistent with the observational 
data can still lead to a recovery of the cosmic ray flux at higher energies than 
presently observed. 


3 LIV Framework 


Coleman and Glashow have proposed a simple formulation for breaking LI 
by a small first order perturbation in the free particle Lagrangian [14]. This 
formalism has the advantages of (1) simplicity, (2) preserving the 51/(3) (8) 
51/(2) (8)17(1) standard model of strong and electroweak interactions, (3) hav- 
ing the perturbative term in the Lagrangian to consist of operators of mass 
dimension 4 that thus preserves power counting renormalizability, and (4) be- 
ing rotationally invariant in a preferred frame that can be taken to be the rest 
frame of the 2.7 K cosmic background radiation. This formalism has proven 
useful in exploring astrophysical data for testing LIV [14], [19], [28]. 

To accomplish this, Coleman and Glashow start with the free particle La- 
grangian 


C = d^*Zd^ - T*M 2 T 


( 1 ) 
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where 'I' is a column vector of n fields with U(l) invariance and the positive 
Hermitian matrices Z and M 2 can be transformed so that Z is the identity 
and M 2 is diagonalized to produce the standard theory of n decoupled free 
fields. 

They then add a leading order perturbative, Lorentz violating term con- 
structed from only spatial derivatives so that 

C -> C + diVe#®, (2) 


where e is a dimensionless Hermitian matrix that commutes with M 2 so that 
the fields remain separable and the resulting single particle energy-momentum 
eigenstates go from eigenstates of M 2 at low energy to eigenstates of e at high 
energies. 

The Lorentz violating perturbative term shifts the poles of the propagators, 
resulting in free particle dispersion relations of the form 

E 2 = p 2 + m 2 + ep 2 . (3) 

These can be put in the standard form for the dispersion relations 

E" = p 2 c 2 MAV + rn 2 c A MAV , (4) 


by shifting the renormalized mass by the small amount m m/( 1 + e) and 

shifting the velocity from c (=1) by the amount Cmav — y^l + e) ~ 1 + e/2. 


Since the group velocity is given by 
d E \p\ 

9 la yJW + mV MAV 


cmav as \p\ 


CX), 


( 5 ) 


Coleman and Glashow thus identify Cmav as the maximum attainable velocity 
of the free particle. 

Using this formalism, it becomes apparent that, in principle, different particles 
can have different maximum attainable velocities (MAVs) resulting from the 
individually distinguishable eigenstates of the e matrix. These various MAVs 
can all be different from c as well as different from each other. Hereafter, we 
denote the MAV of a particle of type i by c t and the difference 

c i - Cj = ^ Sij ( 6 ) 
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There are other popular formalisms that are inspired by quantum gravity 
models or by speculations on the nature of space-time at the Planck scale, 
1/Mpi ~ 1.5 x 10 -35 m, where M Pt = 1 / \J~G ~ 1.2 x 10 19 GeV. Such for- 
malisms, in the context of effective field theory, can be expressed by postu- 
lating Lagrangians containing operators of dimension> 5 with suppression 
factors as multiples of Adpi [19], [29]. This leads to dispersion relations having 
a series of smaller and smaller terms proportional to p n+2 /Mp l ~ E n+2 /Mp t , 
with n > 1. However, in relating LIV to the observational data on UHECRs, 
we find it useful to use the simpler formalism of Coleman and Glashow. Given 
the limited energy range of the UHECR data relevant to the GZK effect, 
this formalism can later be related to possible Planck scale phenomena and 
quantum gravity models of various sorts. 

We now consider the photomeson production process near threshold where 
single pion production dominates, 

p + 7 — > IV + 7T. (7) 


Using the normal Lorentz invariant kinematics, the energy threshold for pho- 
tomeson interactions of UHECR protons of initial laboratory energy E with 
low energy photons of the GBR with laboratory energy u is determined by 
the relativistic invariance of the square of the total four-momentum of the 
proton-photon system. This relation, together with the threshold inelasticity 
relation E n = [m/(M + m)\E for single pion production, yields the threshold 
conditions for head on collisions in the laboratory frame. In terms of the pion 
energy for single pion production at threshold 


4uE n 


m 2 (2M + m) 
M + m 


( 8 ) 


where M is the rest mass of the proton and m is the rest mass of the pion [3] . 

If LI is broken so that c n > c p , it follows from equations (3), (6) and (8) 
that the threshold energy for photomeson production is altered because the 
square of the four-momentum is shifted from its LI form so that the threshold 
condition becomes 


4a )E n 


m 2 (2M + m) 
M + m 


+ 2 S np E 2 


(9) 


Equation (9) is a quadratic equation with real roots only under the condition 

2a; 2 (M + m) 


$irp — 


m 2 (2M + m) 


~ or jm 2 . 


( 10 ) 
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Defining u> 0 = kT = 2.35 x 10 4 eV, equation (10) can be rewritten 


S np < 3.23 x 10” 24 (cu/cc;o) 2 . 


( 11 ) 


If LIV occurs and 8 np > 0, photomeson production can only take place for 
interactions of CBR photons with energies large enough to satisfy equation 
(11). Single photon photomeson production takes is dominated by the A res- 
onance and takes place close to the interaction threshold. This fact, together 
with equation (9) implies that under some conditions photomeson interac- 
tions leading to GZK suppression can occur for “lower energy” UHE protons 
interacting with relatively higher energy CBR photons on the Wien tail of the 
Planck spectrum, but such interactions for higher energy protons, which would 
normally interact with photons having smaller values of u ;, will be forbidden. 
Thus, the observed UHECR spectrum may exhibit the characteristics of GZK 
suppression near the normal GZK threshold, but the UHECR spectrum can 
“recover” at higher energies owing to the possibility that photomeson interac- 
tions at higher proton energies may be forbidden. 


4 Kinematics 


We now consider a detailed quantitative treatment of this possibility, viz., GZK 
coexisting with LIV. We first give the kinematical relations needed to perform 
our calculations in the presence of a small violation of Lorentz invariance. We 
will denote quantities in the proton rest frame by a prime and quantities in 
the cms system of the proton-photon collision by an asterisk. Quantities in 
the laboratory frame are left unprimed. Following equations (3) and (6), we 
denote 


E 2 = p 2 + 2 S a p 2 + m 2 


( 12 ) 


where 5 a is the difference between the MAV for the particle a and the speed 
of light in the low momentum limit (c = 1). 

The cms energy of particle a is then given by 

\fs~a = \J- E 2 - p 2 = S a p 2 +m 2 a (13) 


where, of course, we must have the condition s a > 0. It is important to note 
that, owing to LIV, in the cms where p = 0 the particle will not generally be 
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at rest because 


d E p 

V = ~dp ^ ~E' 


( 14 ) 


We follow Ref. [18] in defining the square of the total rest energy in the cms 
by 


s — E tot Ptot 


(15) 


Then denoting s p to be the square of the initial proton energy in the system 
where the initial proton momentum is zero, it follows that 

s = 2^/s^e + s p (16) 


where we now define e as the energy of the photon in this system. 

Now let us obtain the expression for the modified inelasticity, K, for the 
photopion producing reaction p + 7 — > N + n. Since the inelasticity is defined 
by the fraction of the total energy carried away by the pion, we can relate the 
energy of the emerging proton and pion to the total energy in the laboratory 
system (essentially the initial energy of the proton) by 


E n = KgE p 

E N = (l-K e )Ep (17) 

where Kg is the inelasticity for a given 9 which is the angle between the 
momentum vectors of the photon and the proton in the laboratory system. In 
order to solve for the inelasticity, we calculate the cms energy of the nucleon 
in two different ways. On one hand, we can use the Lorentz transformation of 
the laboratory nucleon energy to relate it to the cms energy: 


E n = 7 *(E n * + /3*p N * cos 9) 

= 7*(En* + f3*\] E n * — s n (E n ) cos 9), (18) 

where the Lorentz factor for the cms frame is the ratio of the total laboratory 
energy E p + u ~ E p to the total cms energy which is given by equation (13) 
and where we now define u to be the observed energy of the CBR photon in 
the laboratory system.. 

On the other hand, we can derive the cms energy of the nucleon from the 
threshold conditions by replacing the masses of the particles with their rest 
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energies as prescribed by equation (13). This yields the relationship 


‘I^/sEn* — S + Sn — Stx 


( 19 ) 


where the quantities sn and can be determined from equation (13) and are 
given by 


Sir — $7r( KqE p ,) 2 + m n 2 

sn = 5jv[(1 - K 0 )E Pi ] 2 + m 2 . (20) 


Here we have replaced p with E since we can exchange momentum for energy, 
given the high Lorentz factor. We can now combine equations (18) and (19) 
to yield a transcendental equation for Kg: 


(1 — Kg)yfs — 


(s + S N (Ko) — S n (Kg)) 


2^s 


+ P\ 


(s + s N (K e )-s n (Ke )) 2 


4 s 


sn cos 6. 


( 21 ) 


The total inelasticity, K, will be an average of Kg with respect to the angle 
between the proton and photon momenta, 6: 


K 


1 

7 r 


7 r 


0 


Kgdd. 


( 22 ) 


The primary effect of LIV on photopion production is a reduction of phase 
space allowed for the interaction. This results from the limits on the allowed 
range of interaction angles implied by equations (21) and (22). As the pion rest 
energy grows, the cosine term in equation (21) becomes larger. For collisions 
with 6 < 7t/2, kinematically allowed solutions become severely restricted. The 
modified inelasticity that results is the key in determining the effects of LIV 
on photopion production. The inelasticity rapidly drops for higher incident 
proton energies. 

As shown in Ref. [14], in order to modify the effect of photopion production 
on the UHECR spectrum above the GZK energy we must have 5 n > S p . It is 
shown in Figure 10 of Ref. [30] that for most of the allowed parameter space 
near threshold 5 n can be as much as an order of magnitude greater than 8 p . 
Therefore, in this paper we will assume that 5 n 3> S p at or near threshold. 
This assumption is also made in Ref. [18]. We will thus take S np ~ S n = 5. 
We have numerically determined that the dependence of our results on the S np 
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Fig. 1. The calculated proton inelasticity modified by LIV for 5 np = 3 x 10 23 as a 
function of CBR photon energy and proton energy. 

parameter dominates over that on the 5 P parameter, as concluded in Ref. [14]. 
The effect of taking a value of 5 P comparable to S n on the UHECR spectrum 
will be presented in a future paper. 

Figure 1 shows the calculated proton inelasticity modified by LIV for a value 
of 5 np = 3 x 1(T 23 as a function of both CBR photon energy and proton energy. 
Other choices for S np yield similar plots. The principal result of changing the 
value of S np is to change the energy at which LIV effects become significant. 
For a choice of 5 np = 3 x 10~ 23 , there is no observable effect from LIV for E p 
less than ~ 2 x 10 20 eV. Above this energy, the inelasticity precipitously drops 
as the LIV term in the pion rest energy approaches m n . 

With this modified inelasticity, the proton energy loss rate by photomeson 
production is given by 

OO 

E* = - ' n|1 - e ~' /2 ™' (23) 


where 77 is the photon threshold energy for the interaction in the cms and er(e) 
is the total y-p cross section with contributions from direct pion production, 
multipion production, and the A resonance. 
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Fig. 2. The calculated proton attenuation lengths as a function proton energy modi- 
fied by LIV for various values of 8 np (solid lines) , shown with the attenuation length 
for pair production unmodified by LIV (dashed lines). From top to bottom, the 
curves are for S np = 1 x 10 — 22 , 3 x 10 — 23 , 2 x 10 -23 ,1 x 1CT 23 ,3 x 10 _24 ,0 (no 
Lorentz violation). 

The corresponding proton attenuation length is given by cE/(dE/dt). This 
attenuation length is plotted in Figure 2 for various values of S np along with 
the unmodified pair production attenuation for comparison. We do not explore 
the effects of modifying pair production through LIV in this paper. 


5 UHECR Spectra with LIV and Comparison with Present Obser- 
vations 


We will start our calculation of LIV modified LIHECR spectra by assuming 
power-law source spectra for the LIHECRs that are chosen to fit the LIHECR 
data below 60 EeV. We then consider the propagation of high energy protons, 
including energy losses resulting from cosmological redshifting, pair production 
and pion production through interactions with CBR photons. 

We shall assume for this calculation a flat ACDM universe with a Hubble 
constant of H 0 = 70 km s -1 Mpc -1 , taking Ha = 0.7 and Q m = 0.3. The 
energy loss owing to redshifting for a ACDM universe is then given by 
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(24) 


-(3 log E/dt)„d, M ft = H 0 \Ja m (l + ' j : + SJ A . 

The attenuation length for protons against pair production is given by 


- (d log E/dt) w = r lp = r n (E) + r e + e -(E), 


(25) 


The attenuation lengths, £ = cE/r(E), for protons against energy loss by both 
pion production, with and without LIV, are shown together with that for pair 
production in Figure 2. The CBR photon number density increases as (l + z) 3 
and the CBR photon energies increase linearly with (1 + z). The corresponding 
energy loss for protons at any redshift z is thus given by 


r JP {E, z) = (1 + z) 3 r[(l + z)E\. (26) 

We take the photomeson loss rate, r n (E), by updating [3] using the latest 
cross sections listed in the Particle Data Group (http://pdg.lbl.gov) and 
in Ref. [31]. We take the pair-production loss rate, r e + e -(E ) from Ref. [32], 

We calculate the initial energy, E^z), at which a proton is created at a red- 
shift z whose observed energy today is E following the methods detailed in 
Refs. [33] and [34], We neglect the effect of possible small intergalactic mag- 
netic fields on the paths of these ultrahigh energy protons and assume that 
they will propagate along straight lines from their source. The total flux of 
emitted particles from a volume element dV = R 3 (z)r 2 drdQ from redshift z 
and distance r with measured energy E is given by 


J(E)dE 


q(Ei, z)dE i dV 
(1 + z)4TrRlr 2 ’ 


(27) 


We assume that the average UHECR volume emissivity is given by q(E i: z) = 
K(z)E~ r . 


We will assume a source evolution q(Ei,z ) oc (1 + z)^ with ( = 3.6, out to a 
maximum redshift of 2.5. This assumption corresponds to a redshift evolution 
that is proportional to the star formation rate. Our results on LIV are insen- 
sitive to the evolution model assumed because evolution does not affect the 
shape of the LIHECR spectrum near the GZK cutoff energy [24,33]. At higher 
energies where the attenuation length may again become large owing to an 
LIV effect, we find the effect of evolution to be less than 10% when compared 
to the no-evolution case (/ = 0). 


Since R 0 — (1 + z)R(z) and R(z)dr = cdt , by integrating equation (27), one 
obtains 
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Fig. 3. Comparison of the HiRes II data with calculated spectra for various 
values of S wp . From top to bottom, the curves give the predicted spectra for 
5^ = 1 x ICR 22 , 3 x 10 — 23 , 2 x ICR 23 , 1 x ICR 23 , 3 x ICR 24 , 0 (no Lorentz violation). 


J(E) 


3cA (0) r 
8v tH 0 


(1 + z)^ 


^/n,„(l + z ) 3 + Ha 



dEi 

dE 


dz. 


(28) 


In this expression, K(0) is determined by fitting our final calculated spectrum 
to the observational UHECR data [24] assuming T = 2.55, which is consistent 
with the Auger data below 100 EeV. 


The results are shown in Figures 3 and 4. 


6 Discussion of Results 


It has been suggested that a small amount of Lorentz invariance violation 
(LIV) could turn off photomeson interactions of ultrahigh energy cosmic rays 
(LIHECRs) with photons of the cosmic background radiation and thereby elim- 
inate the resulting sharp steepening in the spectrum of the highest energy CRs 
predicted by Greisen Zatsepin and Kuzmin (GZK). Recent measurements of 
the UHECR spectrum reported by the HiRes [11] and Auger [12] collabora- 
tions, however, indicate the possible presence of a GZK effect. 

In order to determine the implications for the search for Lorentz invariance 
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Fig. 4. Comparison of the Auger data with calculated spectra for various val- 
ues of Sv p . From top to bottom, the curves give the predicted spectra for 
5 np = lx 10 -22 , 6 x 10 -23 , 4.5 x 10" 23 , 3 x 10" 23 , 2 x 1CT 23 , 1 x 10" 23 , 3 x 1CT 24 , 0 
(no Lorentz violation). 

violation at ultrahigh energies from the analysis of the air shower events ob- 
served by HiRes and AGASA, we undertook a detailed analysis of the spectral 
features produced by modifications of the kinematical relationships caused by 
LIV at ultrahigh energies. In our analysis, we calculate modified UHECR 
spectra for various values of the Coleman- Glashow parameter, S np , defined as 
the difference between the maximum attainable velocities of the pion and the 
proton produced by LIV. We then compare our results with the experimental 
LIHECR data and thereby place limits on the amount of LIV as defined by 
the 5 np parameter. 

Our results show that the amount of presently observed GZK suppression in 
the LIHECR data is consistent with the possible existence of a small amount 
of LIV. In order to quantify this, we determined the value of S np that results 
in the smallest \ 2 for the modeled UHECR spectral fit using the observational 
data above the GZK energy. We find this value to be 4.5 x 10” 23 . We then 
determined the range of acceptable values for S np . This was done by com- 
puting the probablity of getting a \ 2 value at least as small as the \ 2 value 
determined from the fit. We rejected S np values outside of the confidence level 
associated with la. We thus obtained a best-fit range of 8 %p = 4 . 5 I 4 J x 10” 23 , 
corresponding to an upper limit on 8 wp of 6 x 10 ~ 23 , as shown in Figure 4. 

The HiRes spectral data (see Figure 3) do not go to high enough energy to 
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quantitatively constrain LIV. We also note that the Auger spectrum, being 
consistent with no obvious pair-production feature, does not constrain LIV for 
the pair-production interaction. 

A small LIV effect can be distinguished from a higher energy component 
produced by so-called top-down models because the latter predict relatively 
large fluxes of LIHE photons and neutrinos as well as a significant diffuse GeV 
background flux that could be searched for by the Fermi 7 -ray space telescope. 
The Pierre Auger Observatory collaboration has provided observational upper 
limits on the LIHE photon flux that have already disfavored top-down models 
[35]. The upper limits form Auger indicate that UHE photons at best make 
up only a small percentage of the total UHE flux. This contradicts predictions 
of top-down models that the flux of LIHE photons should be larger than that 
of LIHE protons (See Ref. [9] for a review). 

As opposed to the predictions of the top-down models, the LIV effect cuts off 
UHE pion production at the higher energies and consequent UHE neutrino 
and photon production from LIHE pion decay. LIV would also not produce a 
GeV photon flux. 

It is also possible that the apparent modified GZK suppression in the data 
may be related to an overdensity of nearby sources related to the local super- 
galactic enhancement [3 ] . 1 More and better data will be required in order to 
resolve this question. An LIV effect can be distinguished from a local source 
enhancement by looking for LIHECRs at energies above ~200 EeV, as can be 
seen from Figures 3 and 4. This is because the small amount of LIV that fits 
the observational LIHECR spectra can lead to a recovery of the cosmic ray 
flux at higher energies than presently observed. Searching for such an effect 
will require obtaining a data set containing a much higher number of UHECR 
air shower events. 

In the future, such an increased number of events may be obtained. The 
Auger collaboration has proposed to build an “Auger North” array that would 
be seven times larger than the present southern hemisphere Auger array 
(http://www.augernorth.org). Further into the future, space-based tele- 
scopes designed to look downward at large areas of the Earth’s atmosphere as 
a sensitive detector system for giant air-showers caused by trans-GZK cosmic 
rays [36]. We look forward to these developments that may have important 
implications for fundamental high energy physics. 


1 A correlation with nearby AGN has been hinted at in the Auger data [37]. How- 
ever, the HiRes group has found no significant correlation [38]. 
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